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The  major  advantage  of  hydrofoil  craft  is  their  ability  to  sustain 
high  speeds  in  spite  of  severe  sea  conditions.  To  maximize  this  advantage, 
a  reliable  method  is  needed  to  predict  lift  forces  on  hydrofoils  running 
under  waves.  The  results  of  the  present  investigation  show  that  the 
measured  lift  response  of  finite  aspect  ratio  hydrofoils  in  regular  head 
seas  IS  linearly  related  to  the  measured  wave  input.  Thus,  spectral  theory 
can  be  used  with  gust  response  operators  as  well  as  wave  height  spectra 
to  predict  the  statistical  character  of  the  unsteady  lift  on  finite  aspect 
ratio  hydrofoils  in  irregular  head  seas. 

Using  the  measured  gust  response  operators  from  regular  head  sea 
tests  for  two  hydrofoils  with  aspect  ratios  2  and  4,  the  computed  lift 
spectra  were  found  to  be  in  good  agreement  with  the  measured  lift  spectra. 
Reliable  theoretical  predictions  of  the  gust  response  operator  would  facili¬ 
tate  the  application  of  spectral  theory  to  hydrofoil  design.  Accordingly, 
measured  and  predicted  values  of  the  gust  response  operator  are  compared 
for  the  two  models  tested.  Predictions  of  the  gust  response  operators  were 
theoretically  determined  using  the  results  of  two  approximate,  unsteady, 
finite  aspect  ratio,  airfoil  theories:  those  by  Lawrence  and  Gerber  and 
by  Reissner  and  Stevens,  The  measured  and  predicted  magnitudes  of  the 
gust  response  operator  are  in  fairly  good  agreement,  but  the  phase  is  not. 
Some  of  the  discrepancy  is  shown  to  be  due  to  the  approximations  involved 
in  using  oscillating  airfoil  theory  to  predict  the  lift  in  a  traveling  gust. 
It  IS  anticipated  that  accurate  theoretical  predictions  of  the  lift  spectra 
for  a  finite  aspect  ratio  hydrofoil  in  irregular  head  seas  can  be  obtained 
using  an  exact  lifting  surface  theory  to  compute  the  gust  response  operator. 
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INTRODUCTION 

The  major  advantage  of  hydrofoil  craft  is  their  ability  to  sustain 
high  speeds  in  spite  of  severe  sea  conditions,  whereas  displacement  craft 
of  similar  overall  proportions  and  under  the  same  sea  conditions  must 
operate  at  considerably  reduced  speeds.  This  speed  advantage  is  achieved 
by  reducing  the  hydrofoil  craft's  motions  by  controlling  the  lift  forces 
on  the  hydrofoil  system.  In  order  to  maximize  this  advantage  a  reliable 
method  is  needed  to  predict  the  lift  forces  on  the  hydrofoils  under  realis¬ 
tic  conditions  --  including,  for  example,  the  effects  of  finite  aspect  ratio 
and  irregular  seas. 

The  presently  available  analytical  technique  for  predicting  the  sta¬ 
tistical  characteristics  of  ship  motions  in  irregular  seas  is  to  sum  the 
response  of  the  craft  to  each  of  the  harmonic  components  which  compositely 
represent  the  irregular  sea.  In  using  this  approach,  one  presumes  that 
the  ship  motion  is  linearly  related  to  the  wave  input  at  each  frequency. 

Then  for  such  a  linear  system,  the  output  motion  spectrum  can  be  computed 
from  the  wave  input  spectrum  and  the  system  response  operator.  This  re¬ 
sponse  operator  can  be  obtained  either  from  experiments  using  sinusoidal 
input  over  a  range  of  frequencies  or  it  can  be  predicted  theoretically  if 
a  reliable  theory  for  harmonic  input  to  the  system  is  available. 

This  spectral  superposition  procedure  has  been  verified  experi¬ 
mentally  for  the  case  of  displacement  vessels  in  head  seas.^  However, 
it  has  not  yet  been  validated  in  the  case  of  unsteady  forces  on  hydrofoils. 
This  investigation  was  undertaken  to  find  whether  or  not  the  lift  response 
of  fnmte  aspect  ratio  hydrofoils  in  head  seas  could  be  represented  as  a 
linear  system  and,  thus,  to  validate  the  application  of  spectral  theory. 

To  this  end,  two  hydrofoil  models  with  aspect  ratios  2  and  4  were  tested 
in  regular  and  irregular  head  seas. 

Reliable  theoretical  predictions  of  the  gust  response  operator, 
which  is  the  lift  response  to  a  regular  sea,  would  facilitate  the  appli¬ 
cation  of  spectral  theory  to  hydrofoil  design.  Accordingly,  predicted 
values  of  the  gust  response  operators  for  the  two  hydrofoil  models  are 
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compared  with  measured  values.  An  approximate  gust  response  theory  is 
developed  using  two  approximate  oscillating  airfoil  theories.  ’  Although 
more  accurate  unsteady  airfoil  theories  have  been  developed  recently,  as 
well  as  arbitrary  Froude  number  hydrofoil  theories,^  the  computer  programs 
needed  for  the  application  of  these  theories  have  not  yet  been  developed 
for  the  UNIVAC  1105  at  Stevens  Institute  of  Technology. 

This  research  was  sponsored  by  the  Office  of  Naval  Research  under 
Contract  Nonr  263(45). 


APPARATUS  DESCRIPTION 


Two  series  of  tests  were  carried  out:  one  with  an  aspect  ratio  2 
model  and  one  with  aspect  ratio  4,  so  that  the  effect  of  aspect  ratio  on 
the  unsteady  lift  could  be  evaluated  experimentally.  Both  models  were 
constructed  of  brass  and  had  rectangular  planforms  with  thickness  distri¬ 
bution  as  given  by  the  NACA  64^-012  cross-sect i on »  i.e.j  12%  thick  and 
zero  camber.  Each  model  had  a  chord  length  of  3  in.  The  models  were  held 
rigidly  at  zero  angle  of  attack  at  a  depth  of  2  chords  below  the  calm 
water  surface  by  a  12  in.  steel  strut  with  the  same  cross  section  as  the 
hydrof  oils. 

The  strut  extended  vertically  from  the  midspan  of  the  hydrofoil 
to  a  lift  balance  above  water.  The  balance  output  was  calibrated  and 
found  to  be  linear  over  the  range  of  interest  (-20  lb).  The  balance  was 
supported  on  a  carriage  in  Davidson  Laboratory  Tank  No.  3.  The  surface 
wave  height  was  measured  by  a  resistance  wave  wire  at  the  midchord  axis 
of  the  hydrofoil  about  10  in.  outboard  of  the  strut. 

A  dynamic  calibration,  carried  out  with  the  aspect  ratio  2  model 
mounted  on  the  carriage  in  the  tank,  indicated  that  the  magnification 
factor  (apparent  amplification  rate  of  the  actual  input  due  to  flexi¬ 
bility  of  the  supporting  structure)  was  less  than  3%  over  the  frequency 
range  of  0  to  12  cps.  In  addition,  a  spectrum  analysis  of  the  output 
noise  was  carried  out  with  the  foil  operating  at  test  speed  in  calm  water. 
Although  the  lowest  significant  frequency  of  the  noise  spectrum  was  42  cps. 
the  amplitude  was  sufficiently  large  that  the  desired  output  signals  were 
obscured.  Therefore,  passive  low-pass  filters  were  added  in  the  wave  wire 
and  lift  circuits  which  effectively  attenuated  the  noise  level.  The 
attenuation  and  phase  shift  of  the  desired  signals  due  to  the  filters  were 
indeed  appreciable  in  the  range  of  frequencies  of  encounter  of  the  tests. 
However,  in  the  data  analysis,  only  the  relative  magnitude  and  phase  of 
the  lift  compared  with  the  wave  were  needed,  so  that  the  effect  of  the 
filters  was  not  accounted  for  in  the  data  reduction.  Calibrations  were 
carried  out  to  ensure  that  both  filters  had  the  same  characteristics. 
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The  speed  of  the  apparatus  was  determined  from  an  electrical  timer 
which  was  started  and  stopped  as  the  carriage  passed  two  points  in  the 
tank  at  a  known  distance  apart.  The  test  speed  was  held  constant  for  all 
tests  at  15  ft  per  sec  (-.05%  variation)  which  corresponds  to  a  Froude 
number  of  5*28  based  on  full  chord. 

The  filtered  outputs  of  the  lift  balance  and  wave  wire  were  recorded 
on  paper  tape  during  the  regular  sea  tests.  During  the  irregular  sea  tests, 
the  same  outputs  were  recorded  on  both  paper  and  magnetic  tape  to  facili¬ 
tate  data  reduction.  The  continuous  magnetic  tape  records  were  converted 
to  digital  records  on  IBM  cards  by  means  of  an  analog  to  digital  convert¬ 
ing  system  at  Davidson  Laboratory. 

Between  each  regular  sea  test,  a  waiting  period  of  from  five  to 
twelve  minutes  was  allowed  for  the  water  to  calm  down.  The  range  of  wave 
lengths  in  the  regular  sea  tests  gave  reduced  frequencies  of  encounter 
from  0.05  to  0.5.  This  range  covered  the  range  of  frequencies  of  signi¬ 
ficant  energy  in  the  irregular  wave  spectra.  Both  irregular  sea  tests 
were  divided  into  several  runs,  each  of  which  was  obtained  in  a  different 
part  of  the  irregular  sea  pattern. 
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DATA  ANALYSIS 
Regular  Seas 

During  the  regular  sea  tests  with  each  model,  the  lift  amplitude  L, 
phase  it,  surface  wave  amplitude  a,  and  frequency  of  encounter  f^,  were  ob¬ 
tained  from  the  recorded  outputs  of  the  lift  balance  and  wave  wire-  The 
phase  was  determined  as  the  phase  lead  of  the  lift  time  history  with  respect 
to  the  wave  time  history  at  the  midchord  axis.  The  lift  coefficient 
dimensionless  wave  amplitude  r,  at  the  free  surface,  and  the  reduced  fre¬ 
quency  of  encounter  k^,  were  calculated  from  the  defining  relations 

r  =  L  (1) 

''I  2 

2pU  bs 

_  a  (2) 

^  -  b 


where  p  -  density  of  water 
b  =  semi  chord  length 
s  =  semi  span  length 
U  =  model  speed 

The  ratio  of  the  lift  coefficient  to  the  dimensionless  wave  amplitude  was 
obtained  from  these  test  results  at  each  frequency.  This  ratio  and  the 
corresponding  phase  defines  the  gust  response  operator  R(kg),  as 

R(k  )cos  o)  t  =  --t  cos  ((i)  t  +  j^) 
e  e  r  e 

where  u)  =  Zrrf  is  the  circular  frequency  of  encounter.  The  measured 
magnituL  and  phase  of  R(k^)  are  plotted  in  Figs.  1  through  4  for  the 
aspect  ratio  2  and  4  models. 

During  the  aspect  ratio  4  tests,  a  series  of  runs  were  made  at 
approximately  the  same  frequency  of  encounter  but  with  different  wave 
amplitudes.  The  magnitude  and  phase  of  R(k^)  were  computed  as  before  and 
these  results  are  presented  in  Table  I.  This  series  of  runs  provided 
additional  evidence  concerning  the  linearity  of  the  lift  response. 
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Irregular  Seas 

The  sea  surface  and  lift  response  in  an  irregular  head  sea  are 
each  assumed  to  be  composed  of  a  linear  superposition  of  a  large  number 
of  sine  waves  with  random  phase  distribution.  However,  the  magnitude  and 
phase  of  each  component  of  the  lift  response  are  assumed  to  have  a  definite 
relationship  to  the  corresponding  wave  component-  In  fact,  if  the  former 
assumption  is  valid,  then  the  relation  between  the  lift  and  wave  height 
components  at  each  frequency  in  an  irregular  sea  must  be  identical  to  that 
in  a  regular  sea.  The  relation  between  the  lift  and  wave  height  at  various 
frequencies  in  a  regular  sea  is  called  the  gust  response  operator.  Using 
spectral  theory  techniques,  the  statistical  character  of  the  irregular  lift 
can  be  computed  from  the  gust  response  operator  as  well  as  the  statistical 
character  of  the  irre-^ular  sea.  To  verify  the  applicability  of  this  sta¬ 
tistical  method  of  superposition  to  the  case  of  unsteady  lift  on  hydrofoils, 
the  computed  lift  spectra  for  the  two  hydrofoil  models  were  compared  with 
the  measured  lift  spectra. 


The  measured  spectra  were  calculated  by  the  method  outlined  in  Ref..2. 
In  this  method,  the  wave  height  spectrum  is  computed  from  the  time 

history  of  the  wave  r(t),  by 


r ( t ) r ( t-T)dt 


-  i  k  T 
e  dr 


J 


(5) 


the  Fourier  Transform  of  the  auto-correlation  function  (in  square  brackets). 

2 

Similarly,  the  measured  lift  spectrum  ^  was  computed  from  the  time 

history  of  the  lift  C^(t)  by 


C|^(  t  )C|^(  t-T)dt 


-ik  T 
e  ^  dr 


(6) 


The  results  obtained  from  eqs.  5  and  6  for  several  runs  were  then  averaged 
arithmetically  at  each  frequency.  Two  runs  were  used  for  the  Aspect  Ratio 
2  model  and  five  runs  for  Aspect  Ratio  4.  These  average  spectra  are  pre¬ 
sented  in  Figs.  5,  6  and  J  . 
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Under  the  assumptions  outlined  above,  a  computed  lift  spectrum 

Fc  (k  )1  ^  can  be  obtained  from  the  wave  height  spectrum  and  the  measured 
L  L  e  J  c 

gust  response  operator  by  the  relation 


R(k  )  ^  F 

e  L  ^ 


(7) 


Figures  6  and  7  show  the  comparison  between  the  measured  lift 
spectra  obtained  from  eg.  6  with  that  computed  by  means  of  eg.  7  for  the 
two  models. 

A  measure  of  the  linearity  of  the  system  as  well  as  the  guality 
of  the  experiment  is  provided  by  the  coherency,  c(k^).  This  guantity  is 
computed  from  the  relation 


T(k^) 


where  'V{k  )  is  the  cross-spectrum  defined  by 
e 


(8) 


“  2-nJ_ 


1  im 

y  -♦  CD 


r(  t  )C|^(  t-T)dt 


-ik  T 
e  ®  dT 


(9) 


which  is  the  Fourier  Transform  of  the  cross-correlation  function.  The 
coherencies  obtained  from  eg.  8  are  presented  in  Figs.  6  and  7  for  the 
two  sets  of  experiments.  The  measured  phase  relation  between  lift  and 
wave  in  the  irregular  sea  can  be  calculated  by  the  relation 


«(k^) 


tan 


“ Im  v(k  )  “ 
_ ^ 

Re  v(k  ) 
e 


(10) 


The  measured  irregular  sea  phase  angles  are  shown  in  Figs.  2  and  4. 

Since  the  energy  of  a  wave  system  is  proportional  to  the  sguare  of 
the  wave  amplitude,  the  area  under  the  wave  height  spectrum  can  be  used 
to  describe  the  average  properties  of  the  sea.  That  is,  if  R  is  defined  as 
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then  the  average  wave  height  is  given  by 

r  =  0,885  /R  (^2) 

and  this  value  is  reported  in  Fig.  5- 


R-982 
-  9  - 


THEORETICAL  ANALYSIS  --  REGULAR  SEAS 

The  hydrofoil  lift  spectrum  in  an  Irregular  sea  can  be  obtained 
from  eq.  7  by  means  of  the  gust  response  operator.  The  theoretical  approach 
for  predicting  the  gust  response  operator  for  hydrofoils  traveling  at  arbi¬ 
trary  Froude  number  has  been  derived^  but  the  necessary  computer  programs 
for  the  UNIVAC  1105  are  being  developed.  However,  several  approximate 
unsteady  lifting  surface  theories  are  already  available  for  the  case  of 
an  oscillating  airfoil.  Although  these  theories  are  derived  for  the  case 
of  Infinite  fluid,  the  results  can  be  used  in  the  present  study  since  the 
hydrofoils  were  two  chord  lengths  below  the  free  surface.  At  this  depth, 
the  effects  on  the  free  surface  should  be  small. 


Two  approximate  unsteady  airfoil  theories  were  used:  those  by 
3  4 

Relssner  and  Stevens  and  by  Lawrence  and  Gerber.  Each  of  these  theories 
treats  the  problem  of  an  oscillating  wing  in  uniform  flow.  However,  the 
approximate  gust  response  operator  can  be  obtained  from  oscillating  airfoil 
theory  in  the  following  manner. 

The  surface  elevation  T|(x,t)  due  to  a  regular  long  crested  wave  of 
frequency  oj  traveling  in  the  positive  x-directlon  (where  x,y  denotes  an 
orthogonal  Cartesian  coordinate  system  with  the  y-axis  vertically  upward) 

Is  given  by 


T|(x,t)  =  br  cos 


2Tr 


■(U+c)t  j  (13) 


where 

r  =  wave  amplitude  in  semi  chords 

c  =  g/m  =  wave  celerity 
2 

X  =  2TTg/a)  =  wave  length 

U  =  uniform  stream  velocity  In 
positive  x-dl recti  on 

g  =  acceleration  due  to  gravity 
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The  reduced  frequency  of  encounter,  ,  is  then  given  by 


In  eq.  13,  a  crest  is  located  at  x  =  0  (midchord)  when  t  =  0,  which  gives 
the  same  phase  reference  as  was  used  in  the  experiments.  With  the  surface 
elevation  given  by  eq.  13.  the  vertical  velocity  distribution  v(x,y,t)  for 
small  values  of  r  will  then  be 


ZIt 

v(x,y,t)  =  bnoe  ^  T"  "  (U+c)t 


In  the  region  -b  <  x  <  b,  y  =  -bd ,  when  b  «  X,  v(x,-bd,t)  can  be  approxi¬ 
mated  by 


v(x,-bd,t)  =  brwe 


-2tr  .  j 
X  r  2Tfx 


^  cos  UJ  t  -  sin  (D  t 
X  e  e  , 


where  the  relation  oo^  =  2it(U+c)A  has  been  introduced. 

For  a  thin  airfoil  with  zero  camber  oscillating  in  heave  h(t), 
positive  down,  and  pitch  a(t)  about  its  midchord,  positive  for  leading 
edge  up,  such  that 

h(  t )  =  h  b  cos  (0  t  * 

O  6 

O'!  t )  =  <y  cos  (u3  t  +  e ) 


the  vertical  velocity  distribution  over  the  foil  is  given  by 

v(x.-bd.t)  =  +  U  1^)  [h(t)  +  xo.(t)] 

=  (-bh  2Trf  -  U^Qsin  e  -  xry^2TTf^cos  e)  sin  lo^t 

+  (Urv  cos  e  -  xa  2Trf^  sin  s)cos  ui  t  (1 

Comparing  the  corresponding  coefficients  of  the  sine  and  cosine  terms  in 
eq.  15  and  17,  it  is  seen  that  the  vertical  velocity  distribution  on  a 
stationary  airfoil  in  a  sinusoidal  gust  can  be  made  nearly  the  same  as 
that  on  an  oscillating  airfoil  in  uniform  flow  if 
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=  Tr/2 


r  2Trf 


2TTf  b  r 
e 


The  lift  response  to  harmonic  heaving  h{t)  and  pitching  about 
midchord  o'(  t )  can  be  found  from 


r 


Trpb\  ^  2s  Tl  +  L' 

e  L  n  b  O' 


(v{  t ) 


where  L  is  the  complex  lift  amplitude,  s  is  the  semi  span  length  and  L^, 

L'  are  the  complex  unsteady  aerodynamic  derivatives  for  heaving  and  pitching 
O' 

oscillations,  respectively,  with  respect  to  midchord.  The  quantities 
and  h(t)  are  divided  into  real  and  imaginary  parts  in  the  form 

h(t)  =  b(h^  +  ih.)  (cos  ou^t  +  i  sin  ) 

Lu  =  L.  +  iL^,. 
h  hr  hi 

Treating  q'(  t )  and  1^  i  n  a  similar  fashion,  then  expanding  all  terms  in 
the  expression  for  L  e  ®  and  taking  the  real  part  of  both  sides  yields 


RelLe  ®  }  =  -  TTpb^o)  ^  2s  f  Fl.  h  -L.  .+L'  -L'.o  -L'.o'.l 

e  IL  hr  r  hi  ryr  oi  r  O'!  ij 

L,  h.+L.  .h  +L '  cy.+L  .o'  1  sin  1 1 

.  hr  1  hi  r  Q'r  1  rvi  rj  e  J 


cos  t 
e 


The  heaving  motion  is  then  represented  by 


Re  i  ( h  +i  h . )  ( costo  t  +i  sinu)t)^  =  h  cos  out-h.  sinout 
l.'ri  e  ejr  e  i  e 


Comparing  this  result  with  eq,  16,  it  is  seen  that 


h  =  h 
r  o 

h.  =  0 
1 
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Similarly,  for  the  pitching  motion  in  eq.  16  with  e  given  by  eq.  18, 


a  -  0 
r 


a.  -  a ^ 
1  o 


The  lift  response  operator  due  to  harmo^^c  heavi ng  and  pitching  can  then  be 
obtained  from  the  expression  for  Re  jbe  ®  j  giving 

^1  .  V  _  RelLe  J 

—  cos  (o)  t  +  cp) - o 

I"  ®  2pU  bsr 


=  -nk' 


—  L’  .  )  cos  OD  t 
r  Ck"!  e 


h  0? 

L  .  +  —  L'  )sinuDt 
V  hi  r  ctr'  e  J 

The  dimensionless  wave  number  K.  and  the  reduced  wave  frequency  k,  are 
def i ned  as 

2Trb  .  ,  _  (J)b 

K  =  T"  k  -  -g- 

The  approximate  gust  response  operator  is  that  obtained  from  the  lift  re¬ 
sponse  operator  for  the  oscillating  airfoil.  The  amplitudes  of  motion  are 
replaced  by  the  corresponding  wave  characterics  of  eq.  l8.and  the  approximate 
gust  response  operator  is  then 


R(k  )  cos  uu  t  -  -TTkk  e 
0  0  “ 


L  (1  +^)  +  KL' 
hr^  O'! 


cos  UJ  t 
e 


'"-b  -  "  *ir>  '-hr] v} 


The  magnitude  and  phase  of  R(k^)  was  calculated  by  the  methodsof  ref.  3 
and  4  and  the  results  are  presented  in  Figs.  1  and  3  for  1  R(l<g)l  in 

Figs.  2  and  4  for  the  phase  of  R(k^)  for  aspect  ratio  2  and  4.  respectively. 

Using  the  values  of  1  R(l<^)!  predicted  by  eq.  20.  theoretical  pre¬ 
dictions  of  the  lift  spectra,  ^  obtained  from  the  relation 
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[c^(k^)]^  =lR(k^)|2[r(k^)]2  (21) 

The  theoretically  predicted  lift  spectra  are  compared  with  the  measured 
and  computed  lift  spectra  in  Figs.  6  and  7. 
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discussion  OF  RESULTS 
Regular  Seas 


The  measured  and  predicted  values  of  the  magnitude  of  the  gust  re¬ 
sponse  operator  C^/r  are  compared  in  Fig.  1  and  3  for  aspect  ratio  2  and  4. 
respectively.  It  is  seen  that  the  value  of  C^/r  predicted  by  eq.  20  using 
the  Lawrence  and  Gerber  results  is  lower  than  that  predicted  by  Reissner 
and  Stevens.  At  aspect  ratio  2  the  discrepancy  is  about  12%  whereas  at 
aspect  ratio  4,  it  is  only  about  7%* 


In  the  case  of  aspect  ratio  2,  the  value  of  Cj^/r  predicted  by 
Reissner  and  Stevens  is  in  good  agreement  with  the  measured  values  at  all 
frequencies,  while  the  Lawrence  and  Gerber  prediction  falls  below  the 
measured  values,  particularly  at  reduced  frequencies  of  encounter  k^, 
above  0.15.  At  aspect  ratio  4,  the  predicted  value  of  C^/r  using  the 
Reissner  and  Stevens  results  falls  above  the  measured  values  when  0.05<k^<0.3 
whereas  for  0.3<k  <0.5  the  predicted  value  becomes  increasingly  smaller 
than  the  measured  one.  The  Lawrence  and  Gerber  result  at  aspect  ratio  4, 
however,  is  in  agreement  with  the  measured  values  of  Zjr  for  0.05<k^<0.25 
but  again  it  falls  below  them  when  k^>0.25. 


The  measured  and  predicted  values  of  the  phase  cp  of  the  gust  response 
operator  are  compared  in  Figs.  2  and  4  for  aspect  ratio  2  and  4,  respectively. 
It  is  seen  that  the  predicted  values  of  cp  by  the  two  theories  are  in  agree¬ 
ment  but  that  both  predictions  are  less  than  the  measured  values  for  both 
aspect  ratios.  This  discrepancy  is  found  to  be  worse  at  the  high  frequencies 
in  both  cases  and  increases  with  increasing  aspect  ratio. 


At  least  part  of  these  discrepancies  results  from  the  approximations 
used  in  deriving  eq.  20.  An  indication  of  the  magnitude  of  error  intro¬ 
duced  is  seen  by  comparing  the  approximate  gust  response  operator  in 
two-dimensional  flow  with  the  exact  value  as  derived  by  Kemp.  The  values 
of  L  and  L'  in  two-dimensional  flow  are  given  in  Ref.  8. 

h  O' 

The  exact  form  of  the  downwash  due  to  the  regular  waves  is  given 
by  Eq.  14.  However,  to  obtain  the  complex  representation  of  this  downwash, 

consi der 
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•  /  2Tr  , 

i(-^  X  -  U)  t) 


\/(x,t )  -  We 

where  W  is  the  complex  amplitude  of  the  traveling  gust.  Separating  each 
factor  into  its  real  and  imaginary  parts,  then  expanding  and  taking  only 
the  real  part  of  both  sides  leads  to 

Re  •fv(x,t)'l  =  W  cos(-^  X  -  cn  t)  -  W.  si  n(4^x  -  o)  t ) 

(,  J  r  A  e  T  e 

Comparing  this  result  with  eq.  14,  it  is  seen  that 

W  =  0 
r 

2Tr 

W.  =  -b(jore~^  ^ 

1 

The  exact  gust  response  has  been  derived  by  Kemp^  in  the  form 

-  i  (D  t  i  (ij  t ,  p  k 

Le  ®  =  2TTpUbWe  ®  ^  ^ 

where  J  (K)  and  Ji(K)  are  Bessel  functions  and  C(k  )  is  the  Theodorsen 
o  1  e 

function.  If  the  quantity  in  brackets,  which  is  the  generalized  Sears 

function,  is  separated  into  its  real  and  imaginary  parts,  P(k^,K)  and 

Q(k  ,K ) ,  respectively,  then  R„(k  ,K )  is  obtained  as  follows: 
e  be 

RejLe  ®  j- 

R„(k  ,K)cos  (ju  t  =  ~ 

G  e  e  o  .^1 i_ 

2pU  b  s  r 


=  2Trke"'^'^  fQ(k  ,K)cos  w  t  +  P(k  ,K)sin  (u  t|  (24) 
L  e  e  e  e  J 


The  exact  gust  response  operator  obtained  from  eq.  24  and  the  approxi¬ 
mate  gust  response  operator  obtained  from  eq.  20,  are  compared  in  Fig,  8 
for  two-dimensional  flow.  There  the  two  results  are  shown  to  be  in  good 
agreement.  The  error  is  found  to  increase  with  k^  to  about  10%  at  k^  =  0.5* 
However,  the  phases  predicted  by  eqs.  20  and  24  do  not  agree.  The  discre¬ 
pancy  again  increases  with  k^  to  about  22°  at  k^  -  0,5.  Thus,  it  is  con¬ 
cluded  that  the  magnitude  of  the  approximate  gust  response  operator  predicted 
by  eq.  20  is  reliable  in  the  reduced  frequency  range  0<k^<0,5  but  the  phase 
i s  not . 
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However,  these  results  are  for  two-dimensional  flow.  Since  the 
three-dimensional  effect  on  the  gust  response  operator  has  not  as  yet  been 
calculated,  no  quantitative  conclusions  can  be  drawn  regarding  the  accuracy 
of  eg.  20  in  three-dimensional  flow.  However,  in  the  predicted  three- 
dimensional  phase  presented  in  Figs-  2  and  4,  there  is  a  possible  source 
of  error  due  to  the  approximations  leading  to  eq.  20- 

In  order  to  show  the  importance  of  finite  aspect  ratio  in  deter¬ 
mining  unsteady  hydrodynamic  lift,  the  value  of  predicted  by  eq.  20 

for  two-dimensional  flow  is  compared  in  Figs-  1  and  3  with  three-dimensional 
flow  results  for  aspect  ratios  2  and  4,  respectively.  In  both  cases,  there 
is  a  large  discrepancy  near  the  peak  of  the  gust  response  operator .  At 
Aspect  Ratio  two,  the  two-dimensional  value  is  about  100%  higher  than  the 
three-dimensional  value,  while  at  aspect  ratio  four,  the  difference  is 
reduced  to  40%.  In  addition,  the  discrepancy  between  two-  and  three- 
dimensional  predictions  is  found  to  decrease  for  each  aspect  ratio  as  the 
reduced  frequency  of  encounter  increases.  These  observations  are  in  agree¬ 
ment  with  the  conclusions  presented  by  Reissner  and  Stevens^  and  are  further 

substantiated  by  the  graphs  of  L.  ,  L  ,  etc.  given  by  Lawrence  and  Gerber. 

n  Or 

In  addition,  the  same  conclusions  have  been  reached  by  Shiori  and  Tsakonas 
in  their  study  of  unsteady  lifting  surface  theory  applied  to  the  Marine 
propel ler. 

Due  to  the  phase  error  shown  in  Fig.  8,  no  comparison  is  made  here 
between  two-  and  three-dimensional  results.  However,  Shiori  and  Tsakonas 
have  shown  that  the  discrepancy  in  phase  has  the  same  behavior  as  that  for 
the  magnitude  of  the  unsteady  load  on  a  propeller  blade,  i.e.,  the  phase 
of  the  three-dimensional  loading  approaches  that  of  the  predicted  one  in 
two-dimensional  flow  as  the  reduced  frequency  increases  for  a  given  aspect 
ratio,  or  as  the  aspect  ratio  increases  at  constant  reduced  frequency.  It 
is  concluded,  therefore,  that  a  two-dimensional  representation  of  the  un¬ 
steady  forces  is  not  adequate  to  predict  the  magnitude  or  phase  of  the 
lift  on  a  finite  aspect  ratio  hydrofoil  in  head  seas  in  the  reduced  fre¬ 
quency  range  tested  (0<k^<0.5)- 
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Irregular  Seas 

The  averaged  wave  height  spectra  for  the  two  series  of  tests  are 
shown  in  Fig.  together  with  the  average  wave  height  obtained  from  eg.  12. 

The  large  values  of  shown  at  low  may  be  attributed  to  variations 

in  the  spray  sheet  on  the  wave  wire-  The  loss  of  coherency  at  low  as 
well  as  the  low  measured  values  of  this  range  of  k^  shown  in 

Figs.  6  and  7  give  supporting  evidence  for  this  presumption. 

A  computed  lift  spectrum  C  obtained  for  each  aspect 

ratio  from  eg.  7-  The  values  of  I  R(k  )  used  there  were  obtained  from  lines 

e 

faired  through  the  measured  values  of  C^/r  in  regular  seas,  which  are  shown 
in  Figs.  1  and  3-  In  Figs.  6  and  7,  the  measured  and  computed  values  of 
^C^(k^)j^  are  in  good  agreement.  In  addition,  the  values  of  coherency  also 
shown  in  Figs.  6  and  7  are  between  .85  and  1.00  in  the  range  of  0. 08 < k^<  0 . 3  7* 
Finally,  the  measured  values  of  shown  in  Table  I  obtained  at  the  same 

over  a  range  of  r  are  nearly  constant.  From  these  three  observations, 
it  is  concluded  that  the  lift  response  is  linearly  related  to  the  wave  height 
and  thus  the  linear  superposition  theory  using  spectral  technigues  is  appli¬ 
cable  to  this  system. 

The  linear  behavior  of  the  system  indicates  that  nonlinear  effects 
are  not  important.  Thus,  all  approximations  of  linearized  lifting  surface 
theory  are  valid.  One  such  approximation  is  that  the  shed  vorticity  is  a 
plane  sheet  extending  downstream  from  the  trailing  edge.  In  reality,  this 
vortex  sheet  is  distorted  by  the  wave  motion  plus  the  self-induced  velocity 
distribution.  Despite  this,  the  present  results  indicate  that  the  effect 
of  the  distortion  of  the  vortex  sheet  is  small. 

Therefore,  using  linearized  lifting  surface  theory,  a  theoretical 
prediction  of  the  lift  spectrum  in  an  irregular  sea  can  be  obtained  by 
using  eg.  21.  This  prediction  was  carried  out  using  both  the  Lawrence  and 
Gerber  and  the  Reissner  and  Stevens  Theories  in  evaluating  the  approximate 
gust  response  operator.  The  results  are  shown  in  Figs.  6  and  7  for  aspect 
ratio  2  and  4,  respectively.  As  can  be  anticipated  from  the  regular  sea 
results  at  aspect  ratio  2,  the  predicted  lift  spectrum  using  the  Lawrence 
and  Gerber  results  is  lower  than  the  measured  lift  spectrum,  particularly 
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at  the  higher  frequenci es--whereas>  the  Reissner  and  Stevens  prediction 
is  in  fairly  good  agreement  with  the  measured  spectrum.  The  irregular  sea 
results  at  aspect  ratio  4  also  lead  to  the  same  conclusions  as  the  regular 
sea  tests.  That  isj  the  Reissner  and  Stevens  prediction  of  the  lift  spect¬ 
rum  is  too  high  at  low  values  of  k  and  low  at  high  values  of  k  whereas 

G  e 

the  Lawrence  and  Gerber  result  is  in  agreement  with  the  measured  lift 

spectrum  at  low  k  but  too  low  at  high  k  .  However,  in  the  derivation 
“  0 

of  the  unsteady  lift  presented  in  refs.  3  and  4,  approximate  forms  of  the 
integral  equation  relating  the  downwash  to  the  unknown  load  distribution 
were  introduced  in  order  to  simplify  the  numerical  solution.  Now,  by  means 
of  high  speed  computers,  the  numerical  solution  of  the  integral  equation 
is  possible  using  the  exact  form  of  the  kernel  function,  including  the 
free  surface  effect. 

The  measured  phase  of  the  lift  in  the  irregular  sea  with  the  aspect 
ratio  2  model  exceeds  that  obtained  from  the  regular  sea  by  an  amount  which 
increases  with  k^  to  about  15  at  k^  =  0.35<  whereas  at  aspect  ratio  4, 
the  regular  sea  and  irregular  sea  results  are  in  agreement.  This  may  be 
a  result  of  the  smaller  amount  of  data  used  in  the  irregular  sea  analysis 
at  aspect  ratio  2. 

In  this  investigation,  the  hydrofoil  was  considered  to  be  restrained 
against  free  motions.  In  dealing  with  the  dynamic  behavior  of  a  rigid 
hydrofoil  craft  in  a  seaway,  all  restraints  must  be  released  and  the  craft 
must  be  treated  as  a  free  body.  Then,  the  motions  of  the  craft  can  be 
related  to  the  irregular  sea  spectrum  as  was  done  previously  for  the  lift. 
Since  the  same  lift-producing  mechanism  would  be  involved  (response  of  a 
lifting  surface  to  a  downwash  distribution)  in  the  case  of  an  unrestrained 
vehicle  as  in  the  case  of  the  restrained  foil  in  a  traveling  gust,  it  is 
anticipated  that  the  same  conclusions  would  result.  In  particular,  it  is 
expected  that  the  motions  of  a  hydrofoil  craft  in  irregular  head  seas  will 
be  a  linear  superposition  of  the  response  to  the  harmonic  components  of 
the  irregular  sea. 
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CONCLUSIONS 

1.  The  lift  response  of  a  finite  aspect  ratio  hydrofoil  in  head  seas 
is  linearly  related  to  the  wave  height  as  indicated  by  the  good 
agreement  observed  between  measured  and  computed  lift  spectra 

in  irregular  head  seas-  Therefore,  linear  superposition  theory 
can  be  used  to  predict  hydrofoil  lift  spectra. 

2.  Measured  and  predicted  values  of  the  magnitude  of  the  gust  re¬ 
sponse  operator  in  three-dimensional  flow  are  in  fairly  good 
3greement.  However ,  a  discrepancy  was  found  in  the  phase. 

3.  The  predicted  magnitudes  of  the  exact  and  approximate  gust  re¬ 
sponse  operators  in  two-dimensional  flow  show  good  agreement  but 
the  phases  do  not.  Thus,  the  error  introduced  by  the  approxi¬ 
mate  gust  response  operator  in  predicting  the  phase  for  three- 
dimensional  flow  may  account  for  the  discrepancy  with  measured 
va 1 ues . 

4.  In  the  ranges  of  reduced  frequency  of  encounter  and  aspect 
ratio  which  are  typical  for  hydrofoil  craft,  two-dimensional 
theory  is  not  adequate  to  predict  the  gust  response  operator 
in  head  seas. 


ACKNOWLEDGMENTS 


The  authors  wish  to  express  their  appreciation  for  the 
guidance  and  assistance  of  Dr-  Stavros  Tsakonas,  head  of  the 
Fluid  Dynamics  Division  of  Davidson  Laboratory.  Also,  helpful 
discussions  with  Mr.  John  Dalzell  during  the  initial  stages 
of  this  study  is  gratefully  acknowledged. 


R-982 
-  21  - 


REFERENCES 


1.  St.  Denis,  M.  and  Pierson,  W- :  "On  the  Motions  of  Ships  in  Confused 
Seas,"  SNAME  Transactions,  1953- 

2.  Dalzell,  J.  and  Yamanouch i ,  Y.:  "Analysis  of  Model  Test  Results  in 
Irregular  Head  Seas  to  Determine  Motion  Amplitudes  and  Phase  Relation¬ 
ships  to  Waves,"  Ship  Behavior  at  Sea,  Second  Summer  Seminar,  SIT, 

1958. 

3.  Reissner,  E.  and  Stevens,  J.:  "Effects  of  Finite  Span  on  the  Airload 
Distributions  for  Oscillating  Wing,  II  -  Methods  of  Calculation  and 
Examples  of  Application,"  NACA  TN  1195»  19^7- 

4.  Lawrence,  H.  and  Gerber,  E.;  "The  Aerodynamic  Forces  on  Low  Aspect- 
Ratio  Wings  Oscillating  in  an  Incompressible  Flow,"  Journal  of  the 
Aeronautical  Sciences,  Vol .  19>  No-  IK  November  1952. 

5.  Tsakonas,  S.  and  Henry,  C.:  "Three  Dimensional  Gust  Problems  in  the 
Presence  of  Free  Surface  for  Arbitrary  Foil  Orientation  and  Froude 
Number,"  DL  Note  No.  705>  September  1963. 

6.  Spens,  P.:  "A  Digital  Recording  System  for  Model  Tests  in  Irregular 
Waves,"  DL  Note  550>  August  1959- 

7.  Kemp,  N.H.:  "On  the  Lift  and  Circulation  of  Airfoils  in  Some  Unsteady- 
Flow  Problems,"  Journal  of  Aeronautical  Sciences,  October,  1952. 

8.  Scanlan,  R.H.  and  Rosenbaum,  R.:  Aircraft  Vibration  and  Flutter, 

The  Macmillan  Co.,  New  York,  1951- 

9.  Shioiri,  J.  and  Tsakonas,  S.:  "Three-Dimensional  Approach  to  the 
Gust  Problem  for  a  Screw  Propeller,"  DL  Report  940,  March  1963,  to 
be  published  in  Journal  of  Ship  Research. 


R-982 
-  22- 


TABLE  I 

VARIATION  OF  GUST  RESPONSE  OPERATOR 
WITH  WAVE  AMPLITUDE 
ASPECT  RATIO  4 


u 

ft/ sec 

2L 

1b 

2a 

i  n 

f 

e 

cps 

^L 

r 

k 

e 

Cp/r 

deg 

14.98 

14.1 

4.25 

3.63 

.133 

1  .44 

.190 

.0928 

267.8 

14.98 

4.59 

1.32 

3.73 

-0433 

.44 

.195 

.0984 

240-7 

14.98 
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